Selective amplification of scars in a chaotic optical fiber 
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In this letter we propose an original mechanism to select scar modes through coherent gain am- 
plification in a multimode D-shaped fiber. More precisely, we numerically demonstrate the selective 
amplification of scar modes by positioning a gain region in the vicinity of the self-focal point of the 
shortest periodic orbit in the transverse motion. 
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Since their discovery in 1984 by Eric Heller [H], scars 
have remained one of the most intriguing features in the 
field of Wave Chaos. They constitute 'anomalous' modes 
of wave cavities where rays exhibit chaotic dynamics. Ac- 
cording to Random Matrix Theory as well as from 
semiclassical analyses [9] , a spatially uniform distribution 
is expected for modes in such chaotic cavities. Scars devi- 
ate significantly from this ergodic behavior by exhibiting 
unexpected enhancement of intensity in the vicinity of 
the least unstable periodic orbits (PO) of the correspond- 
ing billiards. Since their numerical birth in a quantum 
context, scars have been experimentally observed with 
microwaves [2], with capillary or Faraday waves in 
guided optics ^] and their technological interest is now 
clearly established [l, 0, 0] • The violation of the ergodic 
hypothesis implied by the existence of scars is justified 
by various theoretical approaches ^10|] but remains math- 
ematically uncertain 

In previous works we showed that a passive optical 
fiber with a chaotic-billiard-shape transverse section con- 
stitutes a powerful experimental system for studying (and 
imaging) wave chaos [l^, and particularly when scars 
are involved [4]. Indeed, at the output of a multimode 
D-shaped chaotic fiber we were able to put to evidence a 
scar mode associated to the shortest PO. Unfortunately, 
diffraction due to the finite aperture of the fiber precludes 
any truly selective excitation of a single mode. For an 
input illumination privileging the lowest order modes, a 
few tens of modes contribute to the optical field, and 
the best that one can do is to slightly enhance one of 
the lowest order scars. For larger orders, the increasing 
number of excited modes hampers the scar enhancement 
effect. Scar modes have also been observed in micro- 
lasers [6|. These devices are more and more present in 
photonics, for which high directionality and low lasing 
threshold are required. Scar modes constitute the cav- 
ity modes with the highest quality factors [7] and priv- 
ileged directions of emission. In such lasing microcavi- 
ties, boundary losses play a dominant role in the mode 
selection mechanism and, for a given shape, only a few 
efficient scar modes are supported 0. A different se- 
lection mechanism may be found in the domain of wave 



propagation in random media. In this context, the am- 
plification of localized modes is observed with random 
lasers (see Ref. |13] for a recent review), where they play 
similar to scar modes. Both types of modes rely on strong 
coherent effects and display common statistical features 
[3]. Numerical investigations recently showed that, in 
the localized regime, the lasing modes in a random laser 
just above threshold are the modes of the passive dis- 
ordered system [l5|. The identification between lasing 
and passive modes even holds when the gain region is 
smaller than the mode size[16]. Inspired by these results, 
we present an original way to achieve a selective ampli- 
fication of scar modes in a highly multimode D-shaped 
optical fiber. More precisely, we demonstrate how scar 
modes can be selectively amplified by positioning a gain 
region in the vicinity of specific points along a short PO 
known to give rise to scar modes. We also illustrate the 
robustness of this amplification with a spatially incoher- 
ent speckle-like input. 

We recently showed that multimode D-shaped fibers 
are perfectly suited to achieve high power optical ampli- 
fication [l7|. In the single- mode doped core of an opti- 
cal amplifier, the guided optical field (called the signal) 
is amplified at the cost of an ancillary field (called the 
pump) which propagates in the multimode D-shaped sil- 
ica cladding [l8[. Usually, the optical index of the Er- 
bium or Ytterbium doped core is higher than the index 
of the surrounding D-shaped silica structure. Therefore, 
the signal is guided along the doped region, and, by guid- 
ing the pump, the multimode external core only plays the 
role of an energy reservoir. By some appropriate modifi- 
cations of the fabrication process, which will be presented 
in a forthcoming publication, it is possible to obtain a 
negligible index mismatch between the small active re- 
gion and the rest of the fiber, at least sufficiently small 
to ensure that no mode is guided inside the doped area. 
Thus, both the signal and the pump propagate in the 
entire section of the multimode D-shaped fiber. Here we 
propose a realistic numerical simulation of such a chaotic 
fiber with a localized Ytterbium-doped region. Ytter- 
bium in a silica matrix may be viewed as a four-level sys- 
tem with two metastable levels, which, in our case, can be 
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treated as an effective two-level system. The wavelength 
is Xs = 1020 nm, for the signal, and Xp = 980 nm, for the 
pump, thus limiting the signal reabsorption. The simu- 
lation is based on the Beam Propagation Method (BPM) 
which we already successfully used for simulating 
optical amplification [l7|. 
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FIG. 1: Transverse cross section of the D-shaped multimode 
fiber and 4 scar modes associated to the shortest periodic 
orbit. 

The optical fiber we consider has the D-shaped cross 
section shown in figure [H Its highly multimode core is 
made of pure silica of refractive index rico — 1.451 and 
has a 2R = 125 jim. diameter. The value of the trun- 
cated radius is Rc = R/2. This core is surrounded by 
a polymer cladding of refractive index rid = 1-41 that 
acts both as an optical cladding (to guide light in the 
core) and as a mechanical cladding (to protect the core). 
By using Dirichlet boundary conditions (metallic fiber), 
we are able to solve numerically the Helmholtz equation 
and we show on Fig. [T]some scar modes associated to the 
2-bounce periodic orbit (2P0). As the losses of our sys- 
tem are weak, guided modes of the metallic and dielectric 
fibers are quantitatively the same for the first hundreds 
of modes. The order p of a given scar gives the corre- 
sponding transverse wavenumber hc'^p^ (eigenvalue of Eq. 
(|3]) below) through the following Bohr-Sommerfeld-like 
quantization relation: 



(1) 



where C is the length of the 2P0 and is the phaseshift 
at the core/cladding interface. The additional 7r/2 phase- 
shift is due to the unique self- focal point of the 2P0. The 
active Ytterbium ions are located in a disk of diameter 
15 /im enclosing this self- focal point (see Fig. [1]), that 
also corresponds to a maximum of intensity of the vast 
majority of the 2P0 scars. 

We denote by z the position along the axis of the 
fiber and by f the position in the transverse plane. We 
showed in [j| that the scalar approximation is legitimate 
in the weak guidance limit. Thus the electromagnetic 
field i/jsir^z) representing the signal obeys, in a passive 



fiber, the scalar three-dimensional Helmholtz stationary 
equation 

(A^ + d,,) V^,(r, z) + n\r)k^^s{r. z) = (2) 

where A^ is the transverse Laplacian and ko = 27r/As is 
the vacuum wavenumber. Using the translational invar i- 
ance of the refractive index, Eq. ([2]) can be reduced to a 
stationary Schrodinger equation 
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(/.(f;/3) = -KV(r;/3) (3) 



where (j) is given by iljs{r^z) = / d/3(/)(r; /3)e^^^, the po- 
tential V(r) = (/?^^ - ri^{f)kl)/2 with (3co = ^co^o and 
the transverse wavenumber = (51^ — (5'^ which takes 
on discrete values associated to guided modes (pn- 
Any guided solution ijjsir^z) of Eq. ([2|) can be decom- 
posed on the basis generated by the bound eigenstates 
^n(0 — ^(^5 f^n) of the Schrodinger equation ([3]): 
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,(/)n(r)exp (ipnz) 
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To each mode can be associated an angle On with re- 
spect to the z-axis defined by sin6>n = ^nl f^co- The 
cut-off angle for guided modes is given by sin 6>max = 
^J\ — {nd/ricoY — 0.34, which corresponds to the max- 
imum value tiraax = \/ Pco ~ l^ci ' Kuowiug this value, 
one can obtain the total number of guided modes at 



wavelength Xs which is approximately 3 500 [12|. Mu- 
tatis mutandis^ we can derive similar expressions for the 
electromagnetic field '^^(r, z) of the pump. 

In the paraxial approximation (see [lH for a thorough 
discussion of this approximation in the present context), 
z plays the role of time in a Schrodinger-like equation. By 
using the standard BPM numerical scheme, the coupled 
evolution of the fields over an infinitesimal pseudo-time 
step may be written 



V^5,p(r, z^dz) = exp [-i(D5,p + Ps,p)^z 
+ ^as,p{r, z)dz 



^l^sA^^z) (5) 



where D = — (l/2y^co)^± is the 'diffraction operator' act- 
ing in the Fourier space and P = — (l/2/5co)(^^(^)^o~/^co) 
is the 'propagating operator' accounting for the guided 
propagation in the real space. From the population 
densities Ni{f^z) and N2{r^z) of the two effective Yt- 
terbium energy levels, one obtains the two amplifica- 
tion factors: ap{f^z) = —o-paNi{f^z) and as{r^z) = 
crsa{r]sN2{f,z) - Ni{f,z)), where cisa (cTpa) is the ab- 
sorption cross section of the signal (pump), the stimu- 
lated emission cross section at the wavelength of the sig- 
nal (Jse is also involved through the ratio r]s = (^sel^^sa 
(for a complete account see Ref. [20|, section 14.4). In 
all the results presented below ols ^ 0.1 m~^, preventing 



any significant gain-guiding effect. This effect has been 
recently proposed as a potentiahy effective mechanism 
for providing single-mode operation of high-power fiber 
lasers [22 1. 

In all the simulations performed, the numerical ini- 
tial condition used for the pump evolution ([5]) corre- 
sponds to an illumination by a focussed laser beam, en- 
suring an optimal coupling with a large number of gener- 
ically ergodic modes. The input powers are 6W for the 
pump and 100 /iW for the signal (corresponding to the 
so-called weak signal regime [HI), and the cross sec- 
tions have the following values: apa = 2.65 x 10~^^m^, 
asa = 5.56 X 10-26 = 6.00 x 10"^^ m^. The iden- 

tification of the guided modes which are present in the 
propagating signal is performed through what we call the 
pseudo-time- frequency spectrum C{Hi]z). This quantity 
is derived from a standard procedure in quantum physics 
adapted to the numerical algorithm we use [2l|. More 



precisely, C{hz; z) is obtained from the Fourier transform 
of the correlation function defined as the overlap of the 
propagating field V^s(^, ^) with the initial condition along 
a propagation length L^: 



C{n',z) = j dz' jj dr^ljl{r,{))^lJs{r,z')e-'^^^^'' 

^ 2 

^ ^An{z)8e{K.n- 1^) (6) 
n 

where, for the sake of simplicity, we have written C^n] z) 
in a form which clearly exhibits smoothed delta peaks 
at the values of the transverse wavenumber associated to 
the eigenmodes of the fiber. 

First, we use a selective illumination condition along 
the 2P0 scar mode of order p = 4 which exhibits a 
good overlap with the gain region. Practically, the input 
signal is a plane wave with the transverse wavenumber 
z^^"^) = 10.94/i?. The left inset in Fig. [2] shows the spec- 
trum at 2; = 0. In spite of the plane wave illumination, 
whose wavenumber is associated to the highest peak, the 
unavoidable diffraction generates extra wavenumbers in 
the spectrum. In a passive fiber this spectrum would be 
unchanged during propagation. In Fig. O we show the 
evolution of C^n^z) computed with = 1.3 m and for 
z varying by steps of 0.65 m. The A^^^^-mode is clearly 
enhanced during the propagation. But, as a consequence 
of the initial diffraction, the scar of order p = 3 is also ex- 
cited and amplified. It is worth noting that the spectral 
content of the signal is deeply modified by the amplifica- 
tion: two scar modes completely dominate the spectrum 
at length Zmax = 39 m where the /^^^^-peak reaches a 
plateau. This length corresponds to the maximum am- 
plification length after which the mode is depleted due 
to the absorption of the signal by the active medium. At 
this optimal length, the domination of the scar of order 
j9 = 4 can be directly visualized by imaging either the 
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FIG. 2: Two-dimensional plot of C{Hi; z) defined in (|6|) for 
Lz — 1.3 m. Left inset: spectrum at 2; = 0. Right Inset: 
near-field (NF) and far-field (FF) intensities at the maximum 
signal amplification length Zmax — 39 m. 

We now consider for the signal a spatially incoherent 
illumination resulting from a random superposition of a 
large number of plane waves. We computed the coupled 
evolutions ([5]) for such a random initial condition with 
transverse wave numbers up to = nmax/^- The corre- 
sponding initial signal spectrum was thus almost uniform 
from to R. As depicted in Fig. [3] (a), the evolution of 
C^n^z) exhibits a drastic mode selection along propa- 
gation. Each noticeable peak is associated to a scarred 
transverse wavenumber. This is emphasized in Fig. [3](b) 
where C{n]Zmax) is plotted for Zmax = 45 m and where 
an averaging over five different incoherent illuminations 
has been performed. Almost all scar modes up to or- 
der p = 12 (/^(i2)^ ^ 27.1) are significantly amplified. 
The first peak corresponds to light propagating into the 
fundamental guided mode. The differential amplification 
of the scar modes can be directly related to the relative 
intensities of the modes on the gain zone. This is illus- 
trated in the inset of Fig. [3](b) where the overlap integral 
between the doped area and the intensity of the 2P0 scar 
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modes in a metallic fiber has been plotted. The modes of 
order p = 1, 3, 4, 6 and 10 have a good overlap with the 
gain region and naturally arise in the averaged pseudo- 
intensity spectrum. On the contrary, modes that present 
a bad overlap as those associated to order p = 2, 5 and 7 
are clearly less amplified. 
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FIG. 3: (a) Evolution of C(k; z) along the fiber for an incoher- 
ent initial illumination, (b) Instantaneous spectrum averaged 
over 5 different incoherent illuminations. The inset shows the 
overlap integral between the intensity of the calculated scar 
modes and the doped area. 

In this letter we addressed the problem of scar mode se- 
lection through coherent gain amplification. More specif- 
ically, in a multimode D-shaped optical fiber we demon- 
strated the selective amplification of scar modes by posi- 
tioning a gain region in the vicinity of the self-focal point 
of the shortest periodic orbit in the transverse motion. 
This selection mechanism is so efficient that it can ex- 
tract a family of scar modes from a spatially incoherent 
input signal. Our results may be viewed as the first step 
toward nonlinear wave chaos in guided optics. Indeed, 



pure scar modes could be optically injected in media pre- 
senting cubic nonlinear it ies giving rise to Kerr, four- wave 
mixing or Brillouin effects. Due to its non- uniform spa- 
tial intensity distribution a scar mode would not be ho- 
mogeneously affected by nonlinear dephasing, giving rise 
to an interesting competition between the orderly scar- 
ring phenomenon and destabilizing third-order nonlinear 
effects. 
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